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Abstract 

In this paper, we stndy information transport in mnltiplex networks com¬ 
prised of two conpled snbnetworks. The npper snbnetwork, called the logical 
layer, employs the shortest paths protocol to determine the logical paths 
for packets transmission, while the lower snbnetwork acts as the physical 
layer, in which packets are delivered by the biased random walk mechanism 
characterized with a parameter a. Through simulation, we obtain the op¬ 
timal a corresponding to the maximum network lifetime and the maximum 
number of the arrival packets. Assortative coupling is better than the ran¬ 
dom coupling and the disassortative coupling, since it achieves much better 
transmission performances. Generally, the more homogeneous the lower sub¬ 
network, the better the transmission performances are, which is opposite for 
the upper subnetwork. Finally, we propose an attack centrality for nodes 
based on the topological information of both subnetworks, and further inves¬ 
tigate the transmission performances under targeted attacks. Our work helps 
to understand the spreading and robustness issues of multiplex networks and 
provides some clues about the designing of more efficient and robust routing 
architectures in communication systems. 

Keywords: Traffic dynamics. Routing protocols. Multiplex networks 


1. Introduction 

Communication infrastructures, such as the Internet, mobile phone net¬ 
works, wireless sensor networks etc., become increasingly interconnected. 
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forming large while still growing interdependent communication networks. 
For example, in our daily life people can chat with each other through mo¬ 
bile phone networks, or they can communicate through the Internet by using 
the communication software, such as email, Facebook, etc. Usually the com¬ 
munication flows will travel across many communication networks, if the 
distances between sources and destinations are large. In the battleheld, the 
soldiers and vehicles are connected by different kinds of communication chan¬ 
nels (eg. wireless or not) for reliable communication. In all these systems, 
nodes are connected by more than one kind of links. Each kind of link forms 
a single network layer, and the network layers are often coupled for specihc 
purposes. Beside the technological systems, multiple links or interactions 
among units also widely appear in social and biological systems [H. The 
multiplexity nature of these systems attracts great attention from research 
community in the past few years, and brings about a new research topic, 
namely multiplex networks |2|, l3| . 

Most recently, many researchers focus on the redehnition of the basic 
structure measures coming from single complex network j^, and investigate 
how the multiplexity of layered network structures affect the dynamical pro- 
cesses such as epidemic spreading lB B 0, percolation jsl, 0, [l^, games 


n|, HQ, synchronization [l^, [HI, etc. The transmission of information 
is a common spreading process on communication networks, which is widely 


studied in single complex networks [1^ ll7|, ll8|, ll9|, [20, l2l|]. The most con¬ 


cerned problem is how to alleviate the traffic congestion in order to enhance 
the network capacity. Arenas et al illustrated the phase transition of 


the traffic from the perspective of critical packet generating rate, which was 
thereafter used as the indicator of the network capacity. Zhao et al the¬ 
oretically discussed the network capacity of several traditional topological 
structures by considering the delivery capacity as well as the betweenness 
of the nodes. Most of the other researchers proposed many hard or soft 
strategies to improve the network capacity [1^. The hard strategies are 
modifying the network structure or optimizing the network resources allo¬ 
cation. For example, Liu et al found that it is effective to improve the 
network performance by deliberately removing a small set of edges among 
core nodes. Similarly, Zhang et al obtained that deleting some edges 
among high-betweenness nodes will improve the network capacity. Gong et 
al [2^ obtained that for a given network and a routing strategy, there is an 
optimal resource allocation scheme corresponding to the maximum network 
capacity which is inversely proportional to the average length of the routing 
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paths. 

The hard strategies are too expensive or not applicable in real situation. 
More attention has been paid to the designing of efficient routing strategies 
16l |. In the representative shortest path protocol, traffic flows are directed 


to pass by the core nodes in order to fastly reach the destinations, but this 
will lead to the traffic congestion. Yan et al proposed an efficient routing 
strategy, which intentionally avoids traffic flows passing by the core nodes 
and redistributes the traffic load from core nodes to the marginal nodes. This 
routing strategy greatly improves the network capacity when the delivery ca¬ 
pacity of nodes is identical. Danila et al ^ employed a heuristic algorithm 
to minimize the maximum node betweenness for improving the network ca¬ 
pacity. Wang et al proposed a dynamic routing strategy, which considers 
both the degrees and loads of the neighbor nodes. There are many other ef- 
hcient routing strategies which consider the waiting time of packets in the 
queue j^, the distance of the destinations from the neighbor nodes j^, 32|, 
and so on. 

Until very recently, some researchers studied the traffic dynamics on mul¬ 
tilayer complex networks. Zhuo et al studied the traffic dynamics on 
two-layer complex networks comprised of a logical layer and a physical layer, 
and found that the physical layer is more critical to the overall network ca¬ 
pacity. Du et al ^1 studied the traffic dynamics on coupled spatial networks 
with different travel cost among different network layers. Tan et al in¬ 
vestigated the effects of interconnections on traffic congestion in two coupled 
scale-free networks generated by the Barabasi-Albert (BA) model j^. They 
found that assortative coupling can alleviate traffic congestion much further 
than disassortative coupling and random coupling when node capacity is al¬ 
located according to node usage probability. Nian et al 0 compared several 
routing strategies on two layer degree-coupled networks. Zhou et al stud¬ 
ied the routing issues on multilayered communication networks comprised of 
a wired network and a wireless network. 

In many communication networks, the nodes have a limited power sup¬ 
ply, and they can not work when the batteries run out. Especially in wireless 
sensor networks, the nodes have limited delivery capacity, power, storage 
capacity, and communication distances, and the most important thing is to 
make the networks survive as long as possible. Thus, in the power-limited 

4ll | is more important than network ca- 
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networks, network lifetime 
pacity. In this paper, we study the information transmission in multiplex 
networks comprised of two coupled subnetworks. We focus on network life- 
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time as well as number of arrival packets from the perspective of coupling 
strategies, routing protocols, network structures, and network attacks. 


2. The network model 


In our multiplex networks, there are two network layers. The two network 
layers have the same number of nodes and the same average node degree. 
The topological structures of the network layers are generated by network 
models including the Erdos and Renyi (ER) model ^ and the static (ST) 
model [i^. The ER model generates the random networks with a Poisson 
degree distribution. The static model generates the scale-free networks with 
a power-law degree distribution. Nodes in different layers are connected 
based on one of the three rules: random coupling, assortative coupling, and 
disassortative coupling. For the random coupling, nodes in different layers 
are randomly connected. For the assortative coupling, large-degree nodes in 
one layer are also connected with the large-degree nodes in the other layer. 
For the disassortative coupling, large-degree nodes in one layer are connected 
with small-degree nodes in the other layer. From another point of view, each 
node of the multiplex network has two substitutes, one is in the upper layer, 
and the other is in the lower layer. The two substitutes are connected with 
an edge. 


3. The traffic model 

The upper layer in our multiplex network is the logical layer. In this 
layer, each time every node generates the packets with rate p. If p is 1.5, 
then a node will generate one packet and another one with probability 0.5. 
The destinations of the packets are randomly selected among the other nodes 
except the source nodes. The routing paths for the packets are calculated by 
the shortest path protocol based on the topological structure of the upper 
layer. After the routing paths are available, the packets will be transferred 
to the counterparts of the source nodes in the lower layer for delivery. 

The lower layer is the physical layer for packets delivery. In the lower 
layer, packets are stored in the queue of nodes. The queue length is infinite. 
Each time a node can deliver at most C packets. For simplification purposes, 
we assume that C is infinite to avoid traffic congestion in packets transmis¬ 
sion. When a node wants to deliver a packet, it will check if the destination 
of the packet is among its neighbor nodes. If yes, then the packet will be 
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delivered to the destination node directly, and then be transferred to the 
counterpart in the upper layer immediately. Otherwise, the packet will be 
routed by using the biased random walks strategy. Assume there is a packet 
at node s, then the probability that the packet is sent to the zth neighbor of 
node s is as follow: 


P . = 

^ SI 



( 1 ) 


where ki is the degree of the zth neighbor of node s. a is a control parameter. 
a = 0 means that the packet will be delivered to all the neighbors with equal 
probability. When a > 0, the packet will be sent to the large-degree nodes 
with large probability. When a < 0, the small-degree nodes have a large 
probability to receive the packet. It needs mentioning that all the nodes in 
the routing path calculated in the upper layer should be orderly visited by 
the random walk in the lower layer. We assume that in the lower layer each 
node has E units of energy at the beginning. A node consumes one unit of 
the energy when delivering a packet. We assume that if any node uses up 
its energy, the multiplex network dies. Then, the network lifetime T and the 
total number of arrival packets Q are two representative properties of the 
packets transmission. 


4. The attack centrality 

For monolayer networks, there are many centrality measures for nodes 
such as node degree, closeness, betweenness, eigenvector centrality, PageRank 
centrality, and these measures can be extended to multiplex networks j^. 
Here, we propose a tunable attack centrality based on the degrees of nodes 
in each layer. The tunable attack centrality of node i is as follows: 


2 ;^ = 


^up,i^low,i ’ 


( 2 ) 


where kup^i is the upper-layer degree of node i, and kiow,i is the lower-layer 
degree of node i. /3 is a free parameter which can adjust the relative impor¬ 
tance of kup,i and kiow,i to Zi. When (3 = 0, the centrality value of a node is 
dependent on only its lower-layer degree, while when (3 = 1, the centrality 
value of a node equals to its upper-layer degree. With this tunable attack 
centrality, we can explore how the packets delivery process affected by the 
network targeted attacks. 
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Figure 1: Illustration of our traffic model. The multiplex network is comprised of a logical 
(upper) layer and a physical (lower) layer. A node has a substitute in each of the two 
layers. The two substitutes of a node are connected with an edge marked with dotted 
line (a-a', b-b', c-c', ...). In the upper layer, the routing path between any two nodes is 
calculated with the shortest path algorithm. If there are many shortest paths between two 
nodes, we randomly choose one of them. For instance, the shortest path between node a 
and node c is a—e—c marked with red (online). In the lower layer, the packets are delivered 
with the biased random walk mechanism (Eq. I). For example, there are three possible 
delivery paths from a' to e' which are a' — d! — e' with probability Pa'd'y a' — b' — d! — e' 
with probability Pa'b'Pb'd', and a' — b' — c' — e' with probability Pa'b'Pb'c'- Note that in 
our definition of the random walk, both Pd'e' and Pc'e' equal 1. When delivering a packet, 
a node consumes one unit of the energy. 
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5. Simulation results 

First, we study the influence of the biased random walks on the packets 
transmission on multiplex networks. Specifically, we generate each layer of 
the multiplex networks with the static model. The number of nodes of each 
layer is = 1000 with average node degree (k) = 6 . The power-law pa¬ 
rameter of each layer is 7 = 2.5. The logical layer and the physical layer of 
the multiplex network are interconnected with the one of the three coupling 
strategies including randomly coupling, assortative coupling and disassorta- 
tive coupling. We adjust the control parameter a of the biased random walks 
in the physical layer, and record the corresponding network properties such 
as the network lifetime T and the total number of arrival packets Q. The 
results are shown in Fig. 2, which are the average of 1000 independent runs. 
We see that both T and Q increase with a first, and then decrease with a. In 
the lower layer, when a > 0 , the packets visit large-degree nodes more often, 
and this leads to a fast energy consumption for large-degree nodes, which 
have a relatively small population in the scale-free network. When a < 0, 
the packets are prone to visit small-degree nodes. Although this alleviates 
the traffic loads of large-degree nodes, it will make the packets travel much 
longer in the lower layer, which means more energy will be used in general. 
The optimal a is negative, where T reaches the peak, as shown in Fig. 2. For 
all the coupling strategies, T and Q are positively correlated, and generally 
the longer the network survives, the more packets the network delivers. As¬ 
sortative coupling is the best since it achieves better performances especially 
for Q than the random coupling and disassortative coupling. Disassortative 
coupling is the worst especially for the measure of Q as shown in Fig. 2 
(f). Based on the shortest path protocol, the large-degree nodes have a large 
probability to be on the routing paths of packets. On the other hand, accord¬ 
ing to the random walks theory [i^, the large-degree nodes are more prone 
to be visited by random walkers compared to the small-degree nodes. Thus, 
for the assortative coupling, the packets generally arrive at the destinations 
more early since the large-degree lower layer nodes they often visit happen 
to be in the routing paths, which is opposite for disassortative coupling. 

Then, we study the influence of the network topological structures on 
the traffic dynamics. We mainly focus on the network heterogeneity which 
are controlled by the power-law parameter 7 , and adjust the parameter 7 of 
each layer. The network size N and the average node degree {k) of each layer 
are fixed to be 1000 and 6 respectively. The random coupling, assortative 
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Figure 2: Lifetime T and number of delivered packets Q vs. parameter a of the biased 
random walks for random coupling, assortative coupling and disassortative coupling. Each 
layer of the multiplex network is generated by the ST model with N = 1000, (fc) = 6, and 
7 = 2.5. p = 0.01. E = 3000. Each data point is the average of 1000 independent runs. 


coupling and disassortative coupling are considered for the interconnection of 
the upper layer and the lower layer. We show the results of T as a function of 
7, which are the average of 1000 independent runs. In Fig. 3, we see that for 
random coupling T decreases and converges with increase of 7 of the upper 
layer, when the lower layer is hxed to be either the ER network or the ST 
network. On the other hand, if we fix the upper layer to be either the ER 
network or the ST network, T will increase with 7 of the lower layer. The 
results of Fig. 3 indicate that for random coupling, to achieve large network 
lifetime, the topological structure of the upper layer should be as heteroge¬ 
neous as possible, while the topological structure of the lower layer should be 
as homogeneous as possible. The reason is that a heterogeneous upper layer 
results in a small average routing path, which means generally the packets 
will stay for a short time in the multiplex network and require a relatively 
small amount of energy for their delivery. However, for the lower layer the 
homogeneous network structure leads to a relatively balanced load distribu¬ 
tion, which increases the network lifetime. The results for the assortative 
coupling (Fig. 4) and disassortative coupling (Fig. 5) are similar to that of 
random coupling. There is only an exception for assortative coupling (Fig. 4 
(a)), when the lower layer is fixed to be the ER network. T slightly increases 
with 7, and reaches the peak value around 7 = 3 , and then decreases with 











Lower layer y 


Lower layer y 


Figure 3: Network lifetime T vs. power-law parameter 7 for the random coupling strategy. 
ST on ER means that the upper layer is the ST network and the lower layer is the ER 
network. ST on ST means that both the upper layer and the lower layer are the ST 
networks. ER on ST represents that the upper layer is the ER network, and the lower 
layer is the ST network. For (a) and (b) the lower layer is fixed, while for (c) and (d) 
the upper layer is fixed. Each layer of the multiplex networks have N = 1000 nodes with 
average node degree (k) = 6. p = 0.01. E = 3000. The power-law parameter of the fixed 
ST networks in (b) and (d) is 7 = 2.5. The results are the average of 1000 independent 
runs. 


7. The reason for this exception is that, for assortative coupling the large- 
degree upper layer nodes also have large-degree counterparts in the in the 
lower layer. When the heterogeneity of the upper layer begins to decrease, 
in the lower ER network loads of large-degree nodes are reduced accordingly, 
which results in a direct increase of the network lifetime. When the network 
structure of the upper layer becomes more homogeneous, the average routing 
path is increased accordingly, which leads to a consumption of more energy 
and a decrease of the network lifetime. 

Finally, we study the traffic dynamics in multiplex networks under net¬ 
work disturbance such as the targeted attacks. The attack centrality of each 
node is calculated with Eq. 2. In the attack process, we rank all the nodes 
based on their centrality values and attack the nodes with the largest cen¬ 
trality. It is worth remarking that when we attack a node, we remove its 
counterparts and the incident edges from both the physical layer and the 
logical layer. In the simulation, we mainly focus on how the transmission 
performances such as T and Q affected by the parameter j3 of the attack cen- 
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Figure 4: Network lifetime T vs. power-law parameter 7 for the assortative coupling 
strategy. The multiplex networks and the parameters are the same as in the Fig. 3. The 
results are the average of 1000 independent runs. 






Figure 5: Network lifetime T vs. power-law parameter 7 for the disassortative coupling 
strategy. The multiplex networks and the parameters are the same as in the Fig. 3. The 
results are the average of 1000 independent runs. 
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trality measure. Both the upper layer and the lower layer are generated by 
the ST model. The network parameters of each layer is iV = 1000, {k) = 6 
and 7 = 2.5. The fraction of nodes we attack is set to be 30%. The pa¬ 
rameter a of the biased random walks is set to be 0. The initial amount of 
energy E for each node in the physical layer is 6000. The simulation results 
are shown in Fig. 6, where each data point is the average of 1000 indepen¬ 
dent runs. We see that the results for different coupling strategies are very 
different. For random coupling, T increases first, and then decreases with /5, 
while Q decreases first, and then increases with f3. Those nodes with large 
upper-layer degrees or large lower-layer degrees are the critical nodes for T 
and Q. These critical nodes deliver much more packets than the other nodes, 
and usually use up their energy much early than the others. Thus, removing 
these critical nodes will generally result in a more balanced use of energy, 
which will increase the network lifetime. On the other hand, removing these 
critical nodes on average increases the delivery path length and the difficulty 
of arriving at the destinations for the packets, and this leads to the decrease 
of arrival packets. There is a good compromise between attacking the nodes 
with large upper-layer degrees and the nodes with large lower-layer degrees 
when f3 is around 0.5, where T and Q reach the extreme values, as shown in 
Fig. 6. For assortative coupling, both T and Q are not sensitive to jS. For 
disassortative coupling, there is a phase transition around (5 = 0.5. For all 
the three coupling methods, there are jumps of T and Q when /3 = 1. 

To explain these results, we measure the variation of node rankings based 
on the attack centrality for different (3. We denote the ranking of a node i 
(which has a counterpart in each of the two layer) by r*, which is a function 
of (3. We dehne that the average change of node rankings is as follows: 


(Ar) 


Ef=i(A(/3)-A(0)) 

N 


( 3 ) 


where N is the number of nodes in the multiplex network which equals the 
number of nodes in each of the two layers. We show the results of (Ar) 
in Fig. 7, where the multiplex networks and the parameters are the same 
as in Fig. 6. Each data point is the average of 1000 independent runs. 
For random coupling, (Ar) increases with P continuously, which means the 
node rankings changes continuonsly with /3. For assortative coupling, (Ar) 
is hxed to zero, which indicates that the node rankings do not change, and 
the attacked nodes are actually the same for different (3. For disassortative 
coupling, there is a phase transition of (Ar) aronnd {3 = 0.5, which indicates 
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p p p 

Figure 6 : Network lifetime T and number of arrival packets Q vs parameter /? of the 
attack centrality measure. In the targeted attacks, 30% nodes with the largest centrality 
are attacked. Each layer of the multiplex network is generated by the ST model with 
N = 1000, (k) = 6 , and 7 = 2.5. p = 0.01. E = 6000. Each data point is the average of 
1000 independent runs. 


that the node rankings are constant but different for /3 < 0.5 and (3 > 0.5. 
For all the three coupling methods, there are jumps of (Ar) at /3 = 1, where 
the attack centrality values are calculated by only the topological information 
of the upper layer. 

6. Conclusions 

Multiplex networks theory is an elegant tool for investigating the com¬ 
munication networks which are essentially interconnected with each other. 
In this paper, we study the traffic dynamics embedded in the multiplex net¬ 
works. In previous research of traffic dynamics, the network capacity of 
complex networks is widely discussed. However, we focus on the network 
lifetime and the number of arrival packets, which are more critical than net¬ 
work capacity in power-limited communication networks such as the wireless 
sensor networks. 

In our model, the multiplex networks consists a logical network layer and 
a physical network layer. In the logical layer, the routing paths of packets 
are computed by the shortest path algorithm. In the physical layer, each 
node is power-limited, and packets are transmitted by the biased random 
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Figure 7: (Ar) vs. /3 for different coupling strategies. The multiplex networks and the re¬ 
lated parameters are the same as in Fig. 6. The results are the average of 1000 independent 
runs. 

walk mechanism controlled by a parameter a. 

We obtain the optimal parameters Oopt leading to the largest network 
lifetime and the largest number of arrival packets respectively. We find that 
assortative coupling of the logical layer and the physical layer is much better 
than the random coupling and the disassortative coupling since it achieves 
much better and more robust network performances than the other two cou¬ 
pling strategies. Generally, the more heterogeneous the logical layer, the 
better the transmission performances are, which is opposite for the physical 
layer. 

Finally, we propose a tunable attack centrality with a control parameter 
jS for nodes in multiplex networks, based on which we further study how the 
transmission performances are influenced by the targeted attacks on nodes. 
Simulation results demonstrate that the changes of the network lifetime and 
the number of arrival packets with parameter (3 are very different for dif¬ 
ferent coupling strategies. Generally, our attack centrality is more viable in 
multiplex networks with random coupling and disassortative coupling. 
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